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ABSTRACT 



The composition and amount of interstellar dust within gamma-ray burst (GRB) host galaxies is of key importance when addressing 
selection effects in the GRB redshift distribution, and when studying the properties of their host galaxies. As well as the implications 
for GRB research, probing the dust within the high-z hosts of GRBs also contributes to our understanding of the conditions of the 
interstellar medium and star-formation in the distant Universe. Nevertheless, the physical properties of dust within GRB host galaxies 
continues to be a highly contended issue. In this paper we explore the mean extinction properties of dust within the host galaxies of a 
sample of 17 GRBs with total host galaxy visual extinction Ay < 1 {{Av)=0.4), covering a redshift range z — 0.7 — 3.1. We find 
the average host extinction curve to have an ultraviolet slope comparable to that of the LMC, but with little evidence of a 2175A dust 
extinction feature as observed along Milky Way and LMC sightlines. We cannot at present rule out the presence of a 2 175 A feature, 
and both the standard SMC and LMC extinction curves also provide good fits to our data. However, we can reject an extinction curve 
that has a UV slope as flat as the mean Milky Way extinction curve, whilst also having a 2 175 A feature as prominent as seen in the 
mean Milky Way extinction curve. This is in contrast to the clear detection of a 2175A bump and the flatter extinction curves of some 
more heavily extinguished GRBs (Ay > 1), which may be indicative of there being a dependence between dust abundance and the 
wavelength dependence of dust extinction, as has been previously speculated. 

Key words, gamma-rays: bursts - gamma-ray: observations - galaxies: ISM - dust, extinction 



1. Introduction 

The effects of interstellar dust through its extinction and emis- 
sion properties permeates into most aspects of astrophysics, ei- 
ther as an indirect measure of the environmental conditions and 
enrichment processes present within a galaxy, or as a hindrance 
to be removed from observations. 

The dust extinctio n along numerous lines -of-sight within 
the Milky Way (e.g. Savage & Mathis, 1979) and the Small 
and Large Magel lanic Clouds (SMC and LMC, respectively; 
iFitzpatricki Il989h have been well studied, both in terms of 
the absolute extinction, A\, at a given wavelength A, and 
in terms of the relative Ab — Ay extinction, or reddening, 
E{B — V). The bulk of this work has been focused on de- 
termining the relative extinction as a function of wave l ength , 
or the extinction law or curve fcardelli et al., 198 91: IPeiL Il992t 
iGordon et all 120031: |Fitzpatrick, i2004t iClavtoni i2004 . which 
holds information on the grain size distribution and compo- 
sition of the extinguishing dust along the observed line-of- 
sight. However, it has thus far only been possible to measure 
the dust extinction curve in extragalactic envir onments along 
lines-of-sight to bright background quasars (e.g. Webster et al., 
1995; Pei etal., 1991; Mu rphv & Liske 2004; Maiolino et al., 
2004; El hson et all [200l Iviadilo et aU j2006; Gallerani et al., 
2010) or gamma-ray bursts (GRBs; e.g. lGalama & Wiierst.i200k 



Savaglio et aL 
Starling et al.L 



■20031: IStratta et al.L 12 004: K annetaLl l2006t 
2007t IWatson et al.L l2006>: ,Schadv et al.L 120071) . 



GRBs have been particularly effective at probing the extinc- 
tion cur ve of high-z galaxi e s thanks to t heir in t rinsically simple 
spectra (Chen et al.', 2006; 'Heng et a l.L l2008t iLi et alt l2008bt 
Perleyetal., 2008a, 2010, 2011; ZafareLaU loij- Other ef- 
forts to measure the dust extinction properties in extragalactic 
environments have been made from stacked spectra of large 
samples of gala xies of similar type (e.g. Calzetti et al., 1994 
IViih et all I200I iNoU & PieriiiiL I2005L INoU et al.L 120071. i200% . 
However, such observations are subject to complex radiative 
transfer effects that reprocess the stellar light through numerous 
episodes of gas and dust absorption, emission and scattering. As 
such, the derived galaxy dust extinction properties are highly de- 
pendent on the geometric distribution of the dust, gas and stars 
assumed. The measured wavelength dependent stellar attenua- 
tion is then referred to as the g alaxy's attenuation curve (e.g. 
Calzetti law; ' Calzetti et al.Lll994l) . which is distinct from its dust 
extinction curve (e.g. CCM extinction law; Cardelli et al., 1989j). 

A great deal of progress has been made in the theo- 
retical modelling of the extinction law of dust produced 
by core collapse supe rno yae ( SN e) (e.g. ,Todini & Ferrar; 



2001; 



Nozawa et al. 



,Bianchi & Schneider 



20071), 



20031: [Schneider et d] 
which are thought 



" Hog 

to have 



2 



P. Schady et al.: The Dust Extinction Curve of Gamma-Ray Burst Hosts 



been a possible dominant source of dust formation at high 
redshift (z > 6), when the contribution to dust formation 
from asymptot i c gian t branch stars was still small. Indeed, 
iGallerani et alJ (1201 Oi) found evidence for differences in the 
dust extinction curves in high redshift (z > 4) quasars relative 
to lower redshift quasar extinction curves, suggestive of a 
change in the dust formation or processing mechanism at 
these high redshifts. SN dust extinction curves have also been 
claimed to have been measured in the hosts of high-z q uasars 
(SDSSJ104845.05+463718.3, z=6.193; iMaiolino et all l2004 
and GRBs (GRB 050904, Stratta et al. 2007, although see Zafar 
et al. 2010; GRB 071025, Perley et al. 2010). 

A further ongoing and poorly understood issue is the origin 
and prevalence of the dust extinction feature observed at 2 175 A 
in the Milky Way extinction law, and at lower significance also 
in the Large Magellanic Cloud ( LMC). Small carbona ceous dust 
grains are thought to play a role (iDraine & Leel[T984h . although 
it is still not clear what can produce the observed variation in 
width and height of this feature without varying its central wave- 
length, which shows very small scatter Observationally, a failure 
to account for this feature can also result in a significant under- 
estimation of the am ount of dust extinction in some sources (e.g. 
IConrov et al.L 12010 ^. A more detailed investigation of the varia- 
tion in dust extinction properties with redshift and with environ- 
mental conditions is therefore important to identify the carrier 
of this broad extinction bump, which may have implications for 
dust formation. 

The vast luminosity, simple intrinsic spectra and fading na- 
ture of GRBs make them unique and highly effective probes of 
the extinction properties of dust in high-z gal axies. This is par- 
ticula rly the case for long GRBs (;> 2.0 sec. IKouveliotou et al.L 
[T99l . which are associated with massive star formation (e.g. 
GRB 980425; Galama et al. 1999, GRB 030329; Hjorth et al. 
2003, Stanek et al. 2003), and thus provide lines-of-sight to star- 
forming regions within the host galaxy. Short GRBs typically 
occur far from the centre of their host galaxies, and their faint 
afterglows provide a less effective probe of their (probably low 
density) surrounding environment. In this paper we therefore fo- 
cus on long GRBs only, and all references to GRBs beyond this 
point correspond to this class of burst. 

In this current era of rapid-response GRB observations, the 
number of well-monitored GRBs are now becoming sufficiently 
large to allow statistically meaningful analysis to be done on 
the characteristic properties of GRBs and their host galaxies. 
Early-time, multi-wavelength observations taken with the X-ray 
Telescope (XRT; Burrows et al., 20 051) and the Ultraviolet and 
Optical Telescope (U VOT; Ro ming et al.L 12005 ) onboard Swift 
(iGehrels et al.L l2004l) . in addition to optical and near-infrared 
(NIR) data from rapid-response ground based telescopes, pro- 
vide broadband, high signal-to-noise spectral energy distribu- 
tions (SEDs) for large numbers of afterglows. Dust extinction 
is typically larger at bluer wavelengths, and GRB afterglows 
are, therefore, increasingly attenuated at optical and ultraviolet 
(UV) wavelengths relative to their NIR afterglow. The relatively 
unattenuated rest-frame NIR and X-ray bands (> 2 keV) there- 
fore allow the intrinsic power-law or broken power-law continua 
to be well-determined, from which the imprint left by dust ex- 
tinction can be measured to a high level of accuracy. The blue 
UVOT filters provide unique coverage of the dust-extinction pro- 
file on GRB afterglows at z < 1.5. At higher redshifts, redder 
wavelength observations from ground-based instruments such 
as the Ga mma-Ray Bur st Optical and Near-Infrared Detector 
(GROND; iGreineret all I2OO8I) become important. This longer 



wavelength coverage is also critical for measuring heavily extin- 
guished GRBs which may not be detected by UVOT. 

Typically it is found that for GRBs with well-observed after- 
glow SEDs, the amount of host galaxy extinction is small (to- 



tal T/-band or visual extinc tion (Ay) < 0-3) (e . g. iKann et al, 
200^ ISchadv et all 1201 Ot IStarUng et all l2007t loreiner et al, 



1201 1 ). and has a NIR to X-ray wavelength dependence sim- 
ilar to that of the featureless mean Small Magellanic Cloud 
(SMC). However, recent dedicated GRB rapid-response pro- 
grammes with NIR wavelength coverage are now revealing a 
sample of GRBs (^ 10%) that have bee n significantly extiri - 
guished ( Ay > 1) by host galaxy dust f Greiner et al.L 1201 ih . 
many of which show evidence for a much flatter. Milky Way-like 
host galaxy extinction curve than is observed in the hosts of only 
moderately extinguished GRBs. There are now four GRBs with 
spectroscopically confirmed detections of the 2 175 A absorption 
bump at the redshift of the GRB host galaxy (GRB 070802; 
Eliasdottir et al. 2009, GRB 080607; Prochaska et al. 2009, 
GRB 080605 and 080805; Zafar et al. 2011), as well as sev- 
eral GRBs with a 2 175 A extinction feature detected in the GRB 
SED (GRB 050802; Gates et al. 2007, GRB 070802; Ki-uhler et 
al. 2008, GRB 080603A; Guidorzi et al. 2011). Several claims 
have also been made for GRB hosts having a 'grey' extinc- 
tion law, where the dust extinctiori is only weakly dependent 
on wavelength (e.g. ISava^fio et 120031: IStratta et al'l 120051: 
Perley et al., 2008a!). The causes for this range in the total host 
galaxy dust extinction, and in the shape of the extinction curves 
of GRB hosts are still poorly understood, and better comprehen- 
sion on this issue requires the systematic analysis of a homoge- 
nous data set of GRB afterglows. 



In this paper we explore the variation in the dust extinction 
properties of star-forming galaxies by performing a simultane- 
ous analysis of a sample of 17 NIR to X-ray GRB afterglow 
SEDs. We wish to study the wavelength dependence of the ex- 
tinction curve within GRB host galaxies rather than the absolute 
amount of GRB host galaxy reddening, and we are thus inter- 
ested in analysing the afterglow SEDs over a broad bandpass. 
Our sample has a large redshift range (z — 0.7 — 3.1), which 
provides coverage of the rest-frame wavelength dependence of 
extinction over several decades. As well as addressing the promi- 
nence of a 2 175 A feature and the presence of grey dust in GRB 
host galaxies, the analysis presented in this paper also inves- 
tigates the uniformity of GRB host extinction laws in general. 
The existence of an extinction law specific to GRB host galaxies 
would not only improve GRB afterglow SED modelling and thus 
the accuracy of host galaxy visual extinction measurements, but 
it would also hold information on the grain properties of the dust 
in the GRB surrounding environment. 

In section |2] we give our sample selection criteria and de- 
scribe our data reduction and afterglow SED analysis. Based 
on our individual afterglow SED analysis, we select an optimal 
sample of GRBs with which to study the dust extinction prop- 
erties of GRB host galaxies. This final sample and the conse- 
quent simultaneous SED analysis is described in section |3] Our 
results are presented in sectionlH and a discussion and the main 
conclusions from our work are provided in sections |5] and |6] 
Throughout the paper temporal and spectral indices, a and (3, 
respectively, are denoted such that F{v^t) oc v^^t^"-, and all 
errors are la unless specified otherwise. 
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Table 1. Table listing the 49 GRBs chosen by our selection criteria to have well-sampled afterglow SEDs. 



GRB 


z 


NH,Gal 


Av,Gal 


Epoch^ 


UV/optical/NIR bandpasses 


Rest-frame Band 


Best-fit Model 






(10^^ cm"2) 


(mag) 


(ks) 




Coverage (A) 




050319 


3.24" 


0.13 


0.03 


7 




920-2090 


bknp/MW 


050525A 


0.606'' 


0.91 


0.29 


20 


^4,5j^4 j4,6j4,7^4,8,9^ & U wl m2 w2 


1000-14540 


pow/SMC 


050730 


3.969'^ 


0.30 


0.16 


10 


{KJIry°v b 


790-4700 


pow/MW 


050802 


1.71'* 


0.19 


0.06 


2 


V b u wl Tn2 itj2 


590-2160 


pow/MW 


050820A 


2.615"^ 


0.44 


0.14 


10 


J^^izIRgY^vbuwl 


620-3710 


bknp/LMC 


050922C 


2.198-' 


0.54 


0.32 


5 


R^^~^^v buwl 


710-2360 


pow/MW 


060206 


4.048^* 


0.09 


0.04 


20 


{KHjy^R'^°-^'^vb 


770-4630 


bknp/MW 


Of\Od 1 R 


1 4.0'^ 


u.oo 


(SO 


5 


( KH ^23,24 r25 p26 I, „ ,,,-| ,,,o 


utu — yjo\j 


hknn/SMr 


UDUjUZrt 


i. J 1 


U. 




J 


V b 11 wX 


yuu— juiu 


pu W/ oiVlV^ 


1 1 

uDUj iz 


Z. i 


n 1 s 

V. LJ 


n ns 


J 


i\ s J JX von 


1Z4U— / JjU 


pOW/ olVlL- 


UDUjZD 


J.Zi 


n sn 


n 9 1 

U.Zl 


D.J 


r32,33 p33-40^, 7 


7jU— Z lUU 


po W/ IVl W 


UoUoU/A 


o.Uoz 


U.z4 


o oo 

u.uy 


o 



r-r23 / ■_„\4i I 
_o \^Tg) V 11 


yjU-4zU0 


pOW/bMC 


UoU / ZV 


U.j4 


U.4D 


O 1 T 

U. 1 / 


/O 


% K V u Wl mz wZ 


1040— 30jU 


bknp/bMC 


060904B 


0.703 


1.13 


0.53 


5 


(KJl) Rv buwl m2 w2 


940-13710 


pow/LMC 


060908 


2.43° 


0.23 


0.09 


6 


H bu wl 


660-2200 


pow/SMC 


UoUV 1 z 


u.yj / 


U.jy 


O 1 A 


i 


V b 11 w\ w2 


oJU— jUzU 


pow/MW 


OolUU/ 


1.262' 


0.18 


0.06 


0.7 


(iH) V u wl mz w2 


710-3o50 


pow/bMC 


UDi iZi 






U. 14 





J It V U W 1- JJIZ WZ 




DKnp/ l^iVlL^ 


061126 


1.159" 


1.02 


0.56 


2 


^52^s,53,54^^^^^^2 


740-4030 


bknp/LMC 


070110 


2.352* 


0.16 


0.04 


10 




1150-2250 


bknp/SMC 


070318 


0.836" 


0.14 


0.05 


6 


V bu wl m2 w2 


870-3190 


pow/SMC 


07081 OA 


2.17" 


0.18 


0.07 


5 


R"vbu 


1220-2380 


pow/MW 


071112C 


0.823"' 


0.74 


0.36 


1.3 


(^j)56^57^57-59^ ^57^ ^ ^2 w2 


880-12810 


pow/SMC 


080210 


2.64r 


0.55 


0.26 


5 


KHJz'i'r'g' 


1060-6370 


bknp/LMC 


0803 19B 


0.937" 


0.11 


0.03 


200 


{irf°buwlm2 w2 


830-4490 


pow/SMC 



NOTE- Columns list GRB redshift, Galactic column density, Nn.Cai^ and visual extinction, Av,Gai, in the line-of-sight to the GRB, the 
corresponding SED epoch, the UV/optical/NIR band passes included in the GRB afterglow SED, the rest-frame coverage of the SED, and the 
best-fit SED model. In this last column we indicate both the best-fit continuum (power-law or broken power-law abbreviated to pow and bknp 
respectively) and host galaxy dust extinction model, where we have used the abbreviation MW to indic ate Milky Way extinc tion curve model. 
° iFynbo et al. (2005a); * Folev et al. (2005); " Starling et al. (2005); Fvnbo et al. (20 ()5b|); " |Ledoux et al] (l2005h; ^ Jakobsson et a^ 



b 17^1^,, «f .,1 /onn^A. c ^^ /onn^A. .^i /onncuK. e 'Ledoux et al I j|2005h" ^ 

ai.' ('200^; 

Ledoux et al. (2006); Thoene et al. (2006); " Fusazza et al. ( 2006); ° Rol et al. (2006); ^ lakobsson et al. (2006a); Jakobsson et al. (2006b) 
Bloom et al. (2006); ' Peiiev et al. (2008a); * Jaunse n et all l l2007al) ; " iJaunsen et al.U2007ta) ; " iThoene et al.i l l2007ll ; "' iJakobsson et all ( 120070 : 
Jakobsson et al. (2008a); " Vreeswiik et al. (2008a); 



^ relative to tim e that BAT triggered on the GRB 

^iHuang etalJ 1 2007); ^ Kamble et al. ( 20071); ^IWo zniak et"an d2005h; "^ICobb & Bailvd d2005l) ; '^i Rosenberg & Garnaviclj (l2005h; '^iFlasher et^ 
( 120051); ' lYan agisawa et al. (2005); *^ Homewood et alj (l2005l);~lMirabal et all (|2005|).~ 'Pandey et al. (2006); " Macomb et al. (20O 
Cenk o et al. (2006); Duria c& Price (2005); '.lakobsson et al.' (l2005tf); 'Andreev & Pozanenko (2005); Henvch et al. (2005) 
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17 



iPiranomonte et al. (2005) 
|200] 



Molinari et al. C20070 ; " 



D'Elia et al 

24 



19 



Alatalo et al, 

25 



(2006); 



20 



Sharapov et al j (120061); 
ICovino et al. (2006), Dai et al. (2007) 
(2006) 



(2006); 

31 



Cobb ( 20q6a); 



Curran et al 

26 



(2005) 

Nvsewander et al, 

3° Cenko (2006a); '^^ Milne (2006); ICobbI (2006b); 

Greco etal. (2006); Khamitov et al. (2006) 



21 



(2007); 

KoDoelman (2006); 

33 



Staneketal. (2007); Wozniaketal 
Cenko et al. (2006), '^^ Heartv et al 

34 Tj .:^f ^1 ,onn/;\. 35 



aD 



Ny.sewander et al. (2009); *^ Ouimbv & Rvkoff (2006); 
Mundell etal. (2007) 

54 



IWiersema et al. (2006); 
(120081); '•'^ Smit l ^ et al . 
iDintiniana et al] ( i2007h : 



43 



Cenko (2006b); ' 



Terra etal. (2006); Balivan et al. (2006); 
Morgan & Dai (2006); Rumvantsev et al] 



Cobb&Bailyn (2006); Skvarc (2006) Antonelli et al. (2006); 



46 



Cobb ( 2006q); 

( 20081); Smit l ^ et al (2 006); '^'^ Williams & Milne (2006); Malesani et aL t2007h ; 



59 



iBurenin et ail J2007h ; '^" iRacusin et all d2008h ; 



200^^Uemura et al. (2006); ^'^ GombqcefaD 



lUemura elal., ( ,2007,) ; ,Yoshida et al., ( ,200' 



2. Data Sample and Reduction 

2. 1. Sample Selection 

To study the extinction curve in the local environment of GRBs, 
we require a sample of bursts with afterglow data over a broad 
bandpass, and in particular, covering rest-frame wavelengths 

2000A, where the wavelength dependence of extinction can vary 
significantly between environments (e.g. SMC, LMC and Milky 
Way extinction curves). To carry out our analysis, we select 
those GRBs from the Swift sample up to the end of 2009 that i) 
are classified as long, ii) have a spectroscopic redshift. Hi) were 
observed with the XRT and the UVOT within one hour of trig- 



gering the Swift Burst Alert Telescope (BAT; iBarthelmv et all 
I2005h . iv) had an afterglow detected by both XRT and UVOT 
with a peak UVOT u-band magnitude v < 19, and finally v) 
were detected in a total of at least four UV/optical/NIR filters, 
either ground based or UVOT, where at least one of these filters 
is at a rest-frame wavelength A < 2000A. A spectroscopic red- 
shift is required in order to be able to model the intrinsic GRB 
afterglow spectrum, and we specify the latter three criteria in 
order to select GRBs that have sufficiently high quality broad- 
band afterglow SEDs to constrain our analysis. The X-ray data 
provides a constraint to the afterglow intrinsic broadband spec- 
tral shape, and a further four UV-NIR datapoints are required to 
measure the optical spectral index, /3, the flux normalisation, the 
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Table 1. Continued 



GRB 


z 


NH,Gal 


A-V,Gat 


Epoch 


UV/optical/NIR bandpasses 


Rest-frame Band 


Best-fit Model 






(10^^ cm-2) 


(mag) 


(ks) 




Coverage (A) 




0803 19C 


1.95" 


0.22 


0.08 


0.5 




1310-2260 


pow/SMC 


080411 


1.03' 


0.58 


0.10 


8 


KHJz'i'r'g' 


1900-11420 


bknp/MW 


0804 13B 


1.10^= 


0.31 


0.11 


0.65 


V b u wl m2 w2 


760-2790 


bknp/MW 


080430 


0.767'* 


0.10 


0.04 


10 




910-5010 


bknp/SMC 


080603B 


2.69*^ 


0.12 


0.04 


6.5 


{KHjfR^'^bu 


1050-6330 


pow/SMC 


080710 


0.845^ 


0.41 


0.23 


10 


{KHJzi'r'fv g''^b u wl m2 w2 


870-12560 


pow/SMC 


080721 


2.591" 


0.69 


0.31 


6 


R^°vbu 


1070-2100 


bknp/MW 


080804 


2.204'' 


0.16 


0.05 


1.5 


KHJz'i'r'v g'b u 


1200-7230 


bknp/MW 


080810 


3.35' 


0.33 


0.09 


5 


I^^R^^'^^v g^^b 


900-2030 


pow/SMC 


0809 16A 


0.689^ 


0.18 


0.06 


6 


V b 11 ujI ??i2 w'2 


950-3470 


bknp/SMC 




1 fiQ9''' 
± .uyz 


U. JO 


n 9 1 

U.Z 1 


An 




jyyj — i iU 


nmx/ZSMP 
pUW/vJiVll^ 


no 1 r\r\Q 


i.yoo 


n 9 1 


n 9Q 
V.Iy 


1 
i 


V b u wl 


9An 1 09n 
/ou— ly /u 


pOW/LiVlL- 


ns 1191 

Vol IZl 


9 c:i 9^ 




n 1 ft 
U. 10 


1 n 

iU 


ixiijzivvgou 


1 1 uu— oouu 


DKnp/ oiVlL- 


no 1 TAQ A 

UolzUiA 


z.Uo 


U. 1 / 


U.Uo 




rl3,14 E.13-16,, „15,167 „, 

I ' K V g ' u w 1 


/4U— zyuu 


Dknp/MW 


no 1 TTl 


O ryryO 
2.1 I 


U.zz 


U.Uo 


U. / 


Krtjztrvgouwl 


6UU-OI3U 


DKnp/MW 


non 1 m 
WUIUZ 


1 C /I CP 


U.41 


U. i3 


1 
i 


KHJzirvgo 


15jU-91UU 


DKnp/MW 


uyu^z^ 




n 1 Q 


U.Uo 


1 
i 


yi It) V u w 1 THZ wz 


^^\AC\ ^i^r\ 

lU^U— J 1 J\J 


DKnp/ iVl w 


090618 


0.54'- 


0.58 


0.27 


2 


V b u wl m2 w2 


1040-3800 


bknp/SMC 


090927 


1.37" 


0.29 


0.10 


9 


V b u wl m2 w2 


680-2470 


pow/MW 


091018 


0.971' 


0.28 


0.09 


30 


KHJz'i'r'v g'b u wl m2 w2 


810-11760 


bknp/SMC 


091020 


1.71" 


0.14 


0.05 


6 


V bu wl m2 w2 


590-2160 


pow/LMC 


091029 


2.752" 


0.11 


0.05 


1 


KHJz'i'r'v g'b u wl 


600-6180 


bknp/LMC 


091127 


0.490^" 


0.28 


0.12 


3.7 


HJz'i'r'wl m2 


1340-11950 


bknp/LMC 


091208B 


1.063'= 


0.49 


0.16 


6 


V b u wl m2 w2 


780-2840 


pow/LMC 



" IWiersema et alJ (l2008h: ' iThoene et al.l il20q8ah : IVreeswiik et alJ (l2008d): ICucchiara & Foxl ( 12008 j): " IC ucchiara & Fo xi (l2008bh; ^ 
LPerley et al.' ('2008tf); " 'Jakobsson et al.' ('2008b'); 'Thoene et al.' (I2008cl); ' 'Prochaska et al.l (!2008l); ^ 'Fynbo et al. (2008); ^ ,Vreeswiik et al.l 
(2008c); ' Cucchiara et al. (2008a); Berger & Rauch (2008); " Kuin et al. (2009); ° Cucchiara et al. ( 2008b); de Uaarte Pos tigo et alj(l2009ah: 
« Chornock et al. (2009); Cenko et al. (2009); " Levan et ali ( 120091) ; ' Ide Ugarte Postigo et alj l l2009b() ; " IXu et ali ( |2009|) ; " IChornock et al.l 

(2009);™Cucchiaraetal. (2009);"= Wiersemaetal. (200^ 

^iLietal. (2008a); ^ Thoene et al. (2008b); Terra etal. (2008); Andreev et al. (2008); ^ Oksanen & Hentunen (2008); ^ Miller etal. (2008!); 
'^iKiotz et al. (200 8a); ^.Klotz et al. (2008b); Krti hler et al. (20 09) ^" St arling et al. ( 2009) ; Yoshida et a l, (2008); ,Okuma et al.. t200a) ; 
iRumvantsev et alj ( I2OO8I) ; llsogai & Kawl ( l2008h" ^' ^Voikovl l2008i) ; ^'^ iMori et 'al"J ( [200t) ; ^'^ iKann et al.l j2dlOl) ; 



total host galaxy visual, or V-band (A ^ 5500A) extinction. Ay, 
along the line-of-sight, and the UV slope of the extinction curve. 
This selection criteria leaves a sample of 49 GRBs out of a total 
of 488 GRBs observed by Swift by the end of 2009, and they are 
listed in Table [1] 

By selection, this sample of GRBs will have well-sampled 
afterglow SEDs from which the intrinsic afterglow spectrum 
can be accurately modelled. However, to study the extinction 
curve of their host galaxies, two further selection criteria are 
required. Firstly, there must be evidence for host galaxy dust 
extinction, and secondly, the GRB afterglow should have an in- 
trinsic power-law spectrum from the NIR through to X-ray to 
avoid degeneracy between the afterglow intrinsic spectral slope 
and the amount of dust extinction. This second level of filtering 
requires that we first model the SEDs of each of the initial 49 
GRBs selected above. 

We describe the data reduction and SED analysis on the pre- 
liminary sample of 49 GRBs in sections 1272)12.41 Based on this 
analysis, in section[3]we present our final sample of host galaxy 
dust extinguished GRBs as well as the SED analysis carried out 
to study the average properties of GRB host galaxy extinction 
curves. 

2.2. UVOT, GROND, and ground based data 

The UVOT contains three optical and three UV filters, which 
cover the wavelength range between 160 nm and 600 nm, and 



a clear white filter th at covers the wav elength range between 
160 nm and 800 nm jPoole et al.L l2008h . Photometry was car- 
ried out on pipeline processed sky images downloaded from 
the Swift dat a centre', following the standard UVOT procedure 
dPoole et all l2008i). We applied the tool UV0T2PHA (vl.6) to 
convert UVOT photometry into spectral files compatible with the 
spectral fitting package XSPEC ( Arnaud, 1996), and used version 
V 104 of the 5w//f/U VOT response matrices. In using the UVOT 
filter response curves, we take the known extended tail, or red 
leak in the uvwl (Xc = 2600) and uvw2 (A„ = 1928) filters 
(e.g. lBrown et aUl2010h directly into account. 

For 10 of the 49 GRBs in our sample, data at a similar 
epoch to the Swift observations were available from GROND, 
a seven channel imager covering the 400-2310 nm wavelength 
range simultaneously (g', r', z', z', J, H, Kg) mounted at 
the 2.2 m MPI/ESO telescope at the ESO/La Silla observa- 
tory in Chile. For these 10 GRBs, the GROND data were re - 
duced using the GROND pipeline dKiipcii Yolda^ et all |2008|) . 
and light curves were constructed using standard IRA F tasks 
(Todv, 1993) similar to the procedure described in Kriihle ret al.l 
(2008). GROND optical and NIR photometry was calibrated us- 
ingSDSS (Abazaiian etal ...2009 ) and 2MASS ( Skrutskie et all 
l2006h field stars, respectively. For those GRBs where GROND 
data were not available, optical or NIR photometry redward of 
the UVOT bandpass and reported in refereed papers and GCNs 
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were included in the SED 50% of the sample). For all our 
data, the calibration systemic error was added in quadrature to 
the photometric errors, and where no systematic error was avail- 
able, an error of 0.3 mag was assumed. 

Spectral files were produced for each filter using the ap- 
propriate responsivity c urves. Cousins R and / responsivity 
curves were taken from|Besse ll (119901 and t he J, H and K - 
band respon s ivity cu rves were taken from Cohen et al. l ([T992ah . 
ICohen et al.l (Il992bh and iBessell et all (Il998l) . respectively. In 
constructing these optical response matrices we have not taken 
into account the detector quantum efficiency (QE) curves nor at- 
mospheric effects, which we expect to have a small impact on 
our very broadband spectral modelling. For any non-GROND 
observations taken in the SDSS u'g'r'i'z' filter-set we have as- 
sumed that the response matricies can be well approximated by 
the SDSS-DR6 system response curves^, which do included the 
detector QEs. We expect that the effects of ignoring the small 
differences between the nominal SDSS system and the vari- 
ous true filter-nCCD QEh- atmospheric extinction curves will be 
very small. GROND responsivity curves were taken from the 
GROND home page^ and include the QE of the detectors. 

2.3. XRTdata 

The XRT data were reduced with the Swift XRTPIPELINE tool 
(vO.12.4), and X-ray light curves and spectra in the 0.3-10 keV 
energy range were extracted from the event data using the soft- 
ware package XSELECT (v2.4). Effective area files were created 
with the XRTMKARF tool (vO.5.6), where exposure maps were 
used to correct for bad columns, and response matrices were 
taken from version vOll of the XRT calibration files for both 
WT and PC mode data. The spectral files were grouped to > 20 
counts per energy channel. 

2.4. SED analysis 

We created an X-ray through to optical/NIR instantaneous SED 
for all 49 se lected GRBs following the method described in 
ISchadv et al.l ( 120 IQ) . The epoch of the instantaneous SED was 
selected to be at a time free of spectral evolution, and was also 
chosen to minimise the interpolation in time of the optical and 
X-ray data. Periods of spectral evolution were identified in the 
X-ray data as those where the light curve hardness ratio"* devi- 
ated significantly (90% confidence limit) from the mean. The 
UV/optical/NIR data were screened for spectral evolution by 
verifying the consistency (at 90% confidence) of the best-fit tem- 
poral decay index between filters. 

We fitted all GRB afterglow SEDs within XSPEC (vl2.5.1), 
first modelling the intrinsic continuum with a power-law and 
then with a broken power-law. In the latter case the spectral 
break was modelled as the cooling freque ncy, fc , such that the 
change in slope was fixed at 0.5 (Sari et al.i. ll998h . and the spec- 
tral break was constrained to lie within the UV to X-ray energy 
range. The physical plausibility of such a change in spect ral in- 
dex is also supported by the analysis of IZafar et alJ (1201 ih . who 
did not constrain the spectral slopes in their SED analysis, but 
found A/3 to be typically ^0.5. We considered a broken power- 
law model to provide a better fit to the afterglow spectral con- 
tinuum when the F-test returned a null hypothesis probability 
smaller than 5%. This relatively high probability threshold was 

^ http://www.sdss.org/dr6/instruments/imager/ 
^ http://www.mpe.mpg.de/ jcg/GROND/ 
http://www.swift.ac.uk/xrt_curves/ 



set in order to minimise the chance of any GRBs with an intrinsic 
broken power-law continuum over our observed bandpass being 
incorrectly classified as having a single spectral component. 

The XSPEC model components phabs and zphabs were used, 
respectively, to fit the Milky Way and host galaxy photoelec- 
tric absorption, and zdust was used to model dust extinction in 
the Milky Way and host galaxy. The two photoelectric model 
components translate the soft X-ray absorption resulting from 
medium weight metals into an equivalent neutral hydrogen col- 
umn density using a given absorbing cloud metallicity. We as- 
sumed this metallicity to be solar, both in the case of the 
Milky Way and GRB ho st galaxy, and we us ed the XSPEC de- 
fault solar abundances (lAnders & Grevess5 [1989) . Although 
these abundances ha ve recently been revised to lower values 
(lAsplund et al.'.'2009'). for typical lines-of-sight within the Milky 
Way, the Anders & Grevesse (1989) abundances remain to be 
good estimates (IWatsonil201 1). Also, although GRB host galaxy 
metallicities are typically measured to be sub-solar, the scat- 
ter is large and for many host galaxies, metallicity measure- 
ments are not available. We therefore chose to assume solar host 
galaxy metallicities in our analysis both for consistency within 
our method, as well as to allow us to compare our results with 
those presented in other published literature, where solar abun- 
dances are typically assumed. We fixed the Galactic reddening, 
E{B — V)GaU and the Galactic neutral hydrogen column density 
Nnnai, to the values from the Galactic maps of Schlegel et al] 
d 19981) and lKalberla et al.l(l2005h . respectively. The Galactic dust 
extinction was mod elled on the mean Milky Way extinction law 
as parameterised by iPeil (1 19921) . which has the total-to-selective 
extinction, i?y = Av/E{B — V) set to the average value of 
3.08 measured in the Milky Way diffuse ISM. 

We modelled the dust extinction within the GRB host ga laxy 
on the mean SMC, LMC and Milky Way extinction curves (Te^, 
[1992), which become increasingly flatter blueward of ~ 2500A, 
and have an increasingly pronounced broad extinction feature 
centred at ^2175A (see Fig. [T]i- These three curves thus pro- 
vide good coverage of the range in extinction c urves obse rved 
in the local Universe. To remain consistent with |Peil(ll992h . we 
fixed the total-to-selective extinction of the mean SMC, LMC 
and Milky Way extinction curves to Ry = 2.93, 3.16 and 3.08, 
respectively. This left only the total host galaxy dust reddening 
along the line-of-sight, E{B — V), free to vary. 

Lyman-series absorption in the 912-1215Arest-frame wave- 
length range from neutral hydrogen in th e intergalactic medium 
(IGM) was modelled using results from (lMadauLlT99l . and the 
uncertainty associated to this was added in quadrature to the pho- 
tometric error of any optical data at rest-frame wavelengths blue- 
ward of Lya. The uncertainty on the IGM absorption stems from 
statistical fluctuations in the number of absorbing clouds along 
the line-of-sight, and it is thus a function of redshift, which in 
our sample reached up to 20%. 



3. The GRB mean extinction curve 

To study the dust extinction curve within GRB host galaxies, 
we wish first to reduce our sample of 49 GRBs down to those 
that have both a significant host galaxy dust extinction signature 
imprinted on their SED and that have a single power-law NIR 
to X-ray spectral component. This latter requirement is applied 
because of the degeneracy between the best-fit optical spectral 
index and the steepness of the extinction curve fitted to the SED. 
Over-predicting the steepness of the UV extinction curve when 
modelling the afterglow SED would result in a best-fit optical 
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Table 2. Table listing our final sample of 17 GRBs with afterglows best-fit by a power-law SED model and host galaxy Ay detected 
at 90% confidence. 



GRB 


z 


Host dust 


Av 
(mag) 


(10^^ cm-2) 


P 


X" (dof) 


Null hyp. 
prob. 


050802 


1.71 


mw 


0.49 ±0.10 


2 ■H6+"-'''^ 


0.72 ± 0.04 


73.4 (73) 


0.466 


060502A 


1.51 


smc 






^- ' ^ — 0.05 


23.9 (29) 


0.736 


060607A 


3.082 


smc 


0.08 ±0.04 




0.52 ± 0.02 


29.3 (31) 


0.555 


060904B 


0.703 


Imc 


0.15 ±0.04 


2.98^«-g 


0.85 ±0.01 


36.5 (29) 


0.158 


060912 


0.937 


mw 


0.46t«ii 


^•^^—0 63 


0.85in^ 


16.0 (21) 


0.770 


061007 


1.262 


smc 


0.44 ±0.01 


6.00 ±0.21 


0.89 ±0.01 


285.5 (289) 


0.548 


070318 


0.836 


smc 


0.50 ±0.04 


o r,7+0.60 

O-"'-0.54 


1.14 ±0.02 


63.8 (59) 


0.312 


0708 lOA 


2.17 


mw 


'-'•"^-0.21 




^•'JJ^-0.08 


18.3 (23) 


0.742 


071112C 


0.823 


smc 


20+°-°^ 


0.82ir.^ 


0.58 ± 0.02 


23.7 (17) 


0.128 


0803 19B 


0.937 


smc 


n ofi+^o^ 

^•"0-0.02 


1 -.0+0.42 
^^■^'^-0.37 


0.72 ±0.01 


15.1 (14) 


0.372 


0803 19C 


1.95 


smc 


Q 71+0.08 
'-'■'^-0.07 


Q '7Q + 1-12 

0- ''5-1.03 


0.58 ± 0.04 


99.7 (98) 


0.432 


080710 


0.845 


smc 


0.03 ±0.01 


0.49t«-^ 


0.87 ±0.01 


39.2 (37) 


0.370 


080928 


1.692 


smc 


0.24 ±0.06 


3.54ll:St 


1 09+°"* 

J^-U»-0.05 


33.8 (30) 


0.288 


081008 


1.968 


Imc 


0.29 ±0.07 


7.60tl:g 
2.20+lg 


1 iq + 0.04 
^•^■J-O-OS 


13.2 (23) 


0.947 


090927 


1.37 


mw 




0.88 ± 0.08 


2.2 (9) 


0.989 


091020 


1.71 


Imc 


0.86t«:S° 


r -19+0.96 


1 09+°"^ 
J^-uy-o.os 


40.0 (35) 


0.256 


091208B 


1.063 


Imc 


0.95t«:- 


7.8iti:« 


"■^'-'-0.07 


20.9 (30) 


0.893 



Columns (2) and (3) give the GRB redshift and the best-fit host galaxy extinction law; either smc, Imc or mw to refer to the SMC, LMC or Milky 
Way extinction model fit, and columns (4) to (6) are the best-fit host galaxy Ay, host galaxy Nh, and GRB energy spectral index, /3. The last two 
columns give the and degrees of freedom (dof), and the null hypothesis probability of the best-fit solution. 



spectral index that is flatter than its true value. Such an over- 
prediction would either introduce an artificial break in the after- 
glow broadband continuum, or move the location of a true break 
to smaller energies. In a similar regard, an under-estimation of 
the extinction curve steepness would result in a steeper best-fit 
spectral index than the true value. However, whereas in the for- 
mer case an erroneous best-fit spectral index may still provide an 
acceptable fit to the SED by shifting the position of the cooling 
break, a steeper spectral index would be inconsistent with the X- 
ray spectral index in the context of our SED synchrotron emis- 
sion model. An under-estimation of the extinction curve steep- 
ness would thus result in a poor fit to the data in our SED anal- 
ysis. By only selecting those GRBs with a single intrinsic NIR 
to X-ray afterglow spectral component, we remove this degen- 
eracy, and are thus left with a sample where we know our SED 
analysis to be accurate. This final level of filtering, by definition, 
removes those GRBs from our sample that have host galaxy ex- 
tinction curves with very different slopes to those that we have 
used in our analysis, and we shall return to this point in sec- 
tion |522] 

Of the 49 GRBs analysed, 25 did not require a spectral break 
between the X-ray and optical wavelength ranges. As a consis- 
tency check on our SED analysis, we also verified that the tem- 
poral decay slope in the X-ray and optical bands was the same 
within la errors at the epoch of the instantaneous SED. 36 of 
the 49 GRBs had a host galaxy Ay measured to be greater than 
zero at the 90% confidence, and 17 of these were best-fit by a 
single spectral component, and thus made our final sample. The 
NIR to X-ray SEDs of these 17 GRBs are plotted in Fig. \AJ\ 
with the best-fit FM90 model overplotted, which we described 
in the next section. The SEDs together with the best-fit models 
of the 32 GRBs that did not enter our final sample are plotted in 
Fig. lA.2l In Table |2] we summarise the best-fit model results for 
each of the 17 GRBs in our prime sample. 

To investigate the mean shape of the extinction curve in the 
host galaxies of our sample of 17 GRBs, we applied a simul- 



taneous fit to all 17 NIR to X-ray broadband SEDs using sim- 
ilar models as those described in section 12.41 However, in ad- 
dition to the SMC, LMC and Milky Way extinction laws previ- 
ously fitted to the afterglow SEDs, we also fitted the more gen- 
eral ex tinction curve model developed by Fitzpatri ck & Massal 
to our data, which we shall refer to as the FM90 model. 
This latter extinction curve model has eight free parameters that 
are able to fit variations in the NIR and UV curvature and in 
the prominence of the 2 175 A feature observed along lines-of- 
sight within local and extragalactic galaxies. The large number 
of degrees of freedom in the FM90 model did not make it a 
feasible model to use in the individual afterglow SED analysis. 
However, the large number of data points in our simultaneous 
fit provide sufficient constraint for the FM90 model to be fitted 
to the data. Additionally, the widespread use of this model (e.g . 
Clavton et al., 2003: Go rdon et al .'. 2003: Elfasd ottir et al.Ll2009[ 
Perley et al.,, 201 1) allows us to compare our best-fit results to 
the dust extinction properties of other local and extragalactic en- 
vironments. 

The eight free parameters of the FM90 model are made up of 
the total galaxy dust reddening, E{B — V), the total-to-selective 
dust extinction, i?y, and six further coefficients. These six co- 
efficients define the linear extinction component underlying the 
UV range (ci and C2), the height (C3), width (7), and central 
wavelength (xc) of the 2 175 A bump, which is modelled with a 
Lorentzian-like Dmde profile (D(x\^,Xr)), and the far-UV cur- 
vature (C4). The full equation is given in iFitzpatrick & Mas"sal 
and here we just provide the unexpanded formulism, for 
simplicity: 

A^ = E{B-V)-K (1) 
[Rv ± ci ± + C3D{\-^) + C4F(A-i)] 

The UV intercept, ci, and the UV slope, C2, are tightly cor- 
related, with the most recent anal ysis showing a linear depen - 
dence of ci = 2.09 - 2.89c2 (IFitzpatrick & Massal l2()07h . 
Furthermore, the central wavelength of the 2 175 A bump does 
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Table 3. Best-fit results from our simultaneous SED analysis. 



Model 


— 9 DQ — 9 Rdr-n 


C2 




C4 


'y 

I 




y2 / do f 


Null hyp. 
prob^ 


FM90 


















All data 

^uv > 121;im 


-2.66 
-1.64 


l-31Il.20 


0-00-0.00 
„ ^ .-i_n 57 

0-14Io.i8 


n /fo+0.14 

0- 43Io.ii 

1- 80Io.64 


0.92 
0.92 


4.59 
4.69 


819.9/846 
813.8/814 


0.734 
0.495 


SMC 


















All data 


yf ■in 1 f\ in 
-4.4/±U.ly 


Z.jj±u.1q 




n OT-L/l /I 7 

-u,Zz±U. Ul 


i.uU 


A A/1 

4. ou 


51 J.U/slZ 


U.4u4 


Xuv > 121nm 














786.5/781 


0.438 


LMC 


















All data 


-2.16±0.36 


L31±0.08 


I.92±0.23 


0A2±0.08 


1.05±0.07 


4.63±0.01 


818.7/812 


0.428 


Xuv > 12lnm 














795.4/781 


0.352 


Milky Way 
















5.5e-2« 


All data 


0.12±0.1I 


0.63±0.04 


3.26±0.I1 


0.41±0.02 


0.96±0.01 


4.60=b0.002 


1333.1/812 


Xuv > 12lnm 














1276.3/781 


1.7e-^'^ 



NOTE- Table lists results from both the analysis that used all available afterglow SE D data, and from the analysi s where data covering the Lyman- 
forest were excluded (i.e. no UV data at wavelengths A < 12 1.5 nm) . The best-fit [Fitzpatrick & Massal ( 1 19901) (FM90) parameterisations to the 
SMC, LMC and Milky Way models from lClavton eT^ ( |2003|) are shown in italic as a comparison to the best-fit FM90 fits to our data. 



not vary significantly. We therefore fix Xc to the mean observed 
value of 4.592 fim^^, and tie ci to C2, as described above, leav- 
ing us with six free parameters. Given the typically weak signa- 
ture of a bump in the extinction law of GRB host galaxies to- 
gether with the low spectral resolution available from our broad- 
band data, the constraint on the remaining two bump parameters 
(i.e. C3 and 7) is small, and we therefore also fixed the width 
of the bump, 7, to th e medi an Milky Way value of 7 = 0.922 
(iFitzpatrick & MassaL 12007!) . This allows for a more realistic 
measurement of the bump strength or an upper limit, provided 
by the parameter C3. 

In its original form the FM90 model is only valid up to wave- 
lengths A < 27OOA (rest-frame). Furthermore, most of the GRBs 
in our sample do not have sufficient data at rest-frame NIR wave- 
lengths to derive the total-to-selective extinction, which is give n 
by Rv = 1.1QE{V-K)/E{B-V) (iMorgan & Nandvl[l982l) . 
We therefore appended the mean LMC extinction curve to the 
FM90 analytic model to cover the extinction at wavelengths 
A > 27OOA, and fixed Ry to the mean LMC value, i?y=3.16 
(|PeiLll992l) . We chose the LMC curve because it has a continuum 
that is mid-way between that observed in the SMC and Milky 
Way. Nevertheless, all well-observed extinction curves in the lo- 
cal Universe are comparable redw ard of 2500A, and the spe- 
cific choice of these appended to the'Fitzpatrick & Massal d 19901) 
model at A > 2 700 A has only a marginal effect on our results. 

In our simultaneous SED analysis, all dust extinction models 
fitted had the strength of the host galaxy reddening, E(B-V), as a 
free parameter for each GRB, but all other dust parameters were 
tied between GRBs. The spectral index of the continuum, /3, and 
the amount of host metal absorption, Nh, were also left free to 
vary for each GRB. 

4. Results 

The simultaneous best-fit extinction curves fitted to our data are 
plotted in Fig.[Tl all normalised at A = 3000A. The shaded grey 
area is the 90% confidence region to the best-fit FM90 model. 
As a visual reference we also plot in Fig. |2] all 17 SEDs from 
our sample together with the best-fit model and residuals from 
each of the four host galaxy extinction curve models fitted. A 
summary of our best-fit results is given in Table [3] and to place 
our best-fit FM90 coefficients in context, in this table we also 
list the FM90 best-fit coefficients to the average SMC, LMC and 



Milky Way extinction laws d Clavton et al.L I2003L and references 
therein). 

We find that the Milky Way mean extinction law cannot ad- 
equately fit the host galaxy extinction properties of our sample 
of GRB afterglows, and is rejected at 99.99% confidence (see 
Table [3]). This is also illustrated by the inconsistency between 
the best-fit FM90 extinction curve and the Milky Way mean ex- 
tinction curve plotted in Fig. [T] where the latter is significantly 
flatter and has a more pronounced 2 175 A bump. In the bottom 
right-hand panel of Fig. |2] the residuals from the Milky Way 
model fitted to the data also show a clear excess in the data-to- 
model ratio at the location of the Milky- Way 2 175 A bump. Both 
the SMC and LMC mean extinction laws, on the other hand, 
provide good fits to the data. The best-fit FM90 UV continuum 
is in closer agreement to the LMC rather than the SMC UV 
extinction slope, although at the location of the 2 175 A bump, 
the best-fit FM90 model is closer to the featureless SMC curve. 
Nevetherless, the 90% confidence region of the FM90 model 
plotted in Fig.[T]shows that the presence of a weak 2175A bump 
cannot be ruled out by our data. 

We used the best-fit FM90 extinction curve from our simul- 
taneous SED analysis to re-fit the SEDs of all 49 GRBs that 
were in our initial sample. In these fits, all FM90 parameters 
were fixed to their best-fit value from our simultaneous analy- 
sis, and only the host galaxy reddening, E{B — V), was left 
free to vary. This host extinction curve model provided an ac- 
ceptable fit to all 49 GRB afterglow SEDs, and compared to the 
SMC, LMC and Milky Way models, improved the fit for 65-70% 
of the GRBs with host galaxy dust extinction measured at 90% 
confidence. For the remaining 30-35%, the goodness of fit wors- 
ened by a median of (Ax^) ~ 1 relative to the SMC fit, but by 
(Ax^) = 5 for the Milky Way fit. By allowing the height of 
the 2 175 A bump to vary, the fraction of GRBs better fit by the 
FM90 model rose to 80%. This indicates that for at least 10- 
15% of the extinguished GRBs in our sample, the host galaxy 
extinction curve has a 2 175 A feature. T his is in good agreement 
with the results from lZafar et al.l (1201 ll) . who studied the extinc- 
tion curve for a sample of 41 GRB afterglows using combined, 
broadband photometric and spectroscopic data. In their analysis, 
they found ^ 85% of their extinguished sample to be best-fit by 
SMC-type extinction, and ^ 10% to be best-fit by an extinction 
curve model with a 2 175 A feature. 
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5. Discussion 

In our analysis on the host galaxy dust extinction imprinted on 
the UV through to NIR afterglows of a prime sample of 17 
GRBs, we find the best-fit model to be consistent with both the 
mean SMC and LMC extinction curve. This is in agreement with 
the results from several previo us studies on the dust extinction 



Wavelength (A) 

2000 1500 



prope r ties of GRB hosts (e . g. IStratta et a 



2006; Starli ng et all l2007l TSdiadv et al 



IJ2004; Kannetal. 

, .- 120071: Hengetal., 

2008; Scha dv et aUl2010l)" although co ntradicts the findings of 



Chen et al l (l2006h and lLi et alJ (l2008bl) . who found GRB hosts 



to have a considerably flatter extinction curve, comp arable to or 
flatter than the Milky Way m ean extinction curve. In lChen et al.l 
('2006) and lLi et aP ( 2008bl) the optical intrinsic spectrum, /3opt, 
was inferred indirectly by using the closure relations defi ned 
by synchrotron em ission theory jZhang & Meszarosl |2004|) . In 
IChen et all (l2006h . /3oel 

was calculated from the optical tempo- 
ral decay index, a„pt . iLi eT al. (2008b), on the other hand, used 
a sample of GRBs with consistent X-ray and optical afterglow 
temporal indices, aopt and ax, and assumed the optical spectral 
index, /3opt, to be the same as the measured X-ray spectral in- 
dex, /3x. However, caution must be taken in such analysis, given 
the large fraction of GRBs that have optical and X-ray afterglow 
data that are inconsistent with the sim plest fireb all scenarios (e. g 
iPanaitescu et"all 120061; lEvans et al.L 12009.; .Nardini etall l201d) . 
and although /3opt = Px ir nplies aont = a^, the opposite does 
not necessarily apply (e.g. 'Filg as et all 1201 lb . Further compli- 
cations arise from the degeneracy that there exists between the 
X-ray spectral index and the amount of soft X-ray absorption. 
For example, an error of 0. 1 on the X-ray photon spectral index. 
Fx, would introduce an uncertainty of ^ 0.5 magnitude in the V- 
band flux, and a break between the X-ray and optical emission 
coiTesponding to the cooling frequency (Fopt = Fx — 0.5) would 
reduce the predicted flux in the optical by a factor of ^ 10. 

The shape of our best-fit GRB host extinction curve strongly 
contrast with the extinction curves of the highly extinguished 
GRBs 070802 and 080607 measured in 'EKasddttir et al.^ (l2009l) 
and Perley et al (2011) resp ectively, an d GRB 080605 and 
GRB 080805 as measured by IZafar et alJ (l201lh . all of which 
are flatter and have a clear 2 175 A feature. None of these GRBs 
were included in our sample due to the faintness of their ob- 
served optical afterglows. To investigate the possible origins for 
these differences in host extinction curves we consider first the 
systematic uncertainties present in our analysis, as well as the 
selection effects in our sample. We then compare properties of 
those GRBs in our sample with those of GRBs with distinct host 
dust extinction properties. 



5.1. Systematic biases 

There are two principle areas where systematics may be intro- 
duced, and this is i) in our modelling of the Galactic dust extinc- 
tion and soft X-ray absorption, and ii) in the effect of host galaxy 
neutral hydrogen on our UV and optical data. 



5.1 .1 . Milky Way dust extinction 

The uncertainty in the NH,Gai and Ay.cai can be up to a fac- 
tor of two when the line-of-sight is close to the Galactic plane. 
However, probably as a consequence of our selection criteria re- 
quiring a GRB peak v-band afterglow v < 19, only three of 
the initial 49 GRBs analysed have a Galactic visual extinction 
Ay > 0.5, at which point the uncertainty in the Galactic redden- 
ing becomes significant. Only one of these three GRBs entered 



FM90Rv=3.16 

--■ FM90;noUVdata<1215A 
-■- SMC 
■-■ LMC 

Milky Way 



< 

<; 




23456789 

Inverse wavelength (l/^m) 

Fig. 1. Best-fit mean GRB host extinction law as derived from 
our simultaneous SED fits normalised at A » and corre- 

3000A 

sponding 90% confidence region (grey). The best-fit FM90 
model is plotted in black (solid), and the best-fit model to the 
SEDs with UV data blueward of 12 15 A omitted is shown in 
blue (long-dashed). Also shown for comparison are the mean 
SMC (red, dot-dash), LMC (green, dot-dot-dash) and Milky Way 

I pink, dott ed) extinction laws, using the parameterisation from 
PdL ll992h . 



our final sample. Furthermore, due to the general inverse wave- 
length dependence of extinction, any extinction in the observed 
optical bands is dominated by the uncertainty in the much larger 
UV rest-frame extinction. Therefore, although the SED analy- 
sis of these three more heavily Milky-Way extinguished GRB 
afterglows may be uncertain, this should only introduce scatter 
rather than alter our results. In the case of the Galactic neutral 
hydrogen column, Nn^oai, the uncertainty is typically less than 
10%, which is small relative to the uncertainty on the best-fit 
host galaxy soft X-ray absorption. 

Another source of uncertainty is in the Milky Way extinc- 
tion la w, which varies along different lines-of-sight as a function 
of i?v (ICardenietall[T989l) . The scatter around the mean value 
of Rv=3.l is, however, small, with a root mean square of 0.1 
(Schultz & Wiemer, 1975). Furthermore, larger values of Ry are 
seen along lines-of-sight that cross dense regions of dust, and are 
therefore not very common. This is especially true for our sam- 
ple, which is biased against heavily extinguished lines-of-sight. 
We therefore do not expect variations in the Galactic extinction 
curve along the lines-of-sight to our sample of GRBs to have a 
significant affect on our final results. 



5.1.2. Host neutral hydrogen absorption 

A factor that needs to be addressed before discussing our results 
is the effect of gas absorption at UV wavelengths within the host 
galaxy, which is unaccounted for in our SED modelling. The 
presence of a strong damped Lyman-a (DLA) system within the 
GRB host galaxy, which can contain column densities of neu- 
tral gas as lar ge as Nh ^ 10^^ cm~^ (e.g. GRB 050401 and 
GRB 060926; iFvnbo et al.L l2009l) . may systematicafly bias our 
dust extinction measurements at rest-frame wavelengths 1150- 
1250A. If indeed significant, then this neutral hydrogen absorp- 
tion within the host galaxy may be incorrectly attributed to dust 
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Fig. 2. UV/optical/NlR afterglow SEDs of the 17 GRBs in our sample plotted in the rest-frame. Overplotted are the best-fit models 
from the simultaneous fits, where each of the four figures correspond to a different host galaxy extinction curve model fitted to the 
data. The four extinction curves fitted are our best-fit FM90 extinction curve (top left; best-fit coefficients given in TableO, the mean 
SMC (top right), mean LMC (bottom left), and mean Milky Way (bottom right) extinction curves, where the latter three extinction 
curves use the parameterisation from Pei ( 1992 ). In all panels the y-axis is arbitrary and the afterglow SEDs are mostly ordered 
by afterglow host galaxy extinction with the most extinguished at the bottom, and the least plotted at the top. In all four panels the 
dashed lines indicate the Lyman-break at 912. 5 A, and the location at 1215A, blueward of which the Lyman-forest presents a further 
absorption component. 
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extinction, and thus result in an apparent steepening of the best- 
fit extinction curve in the rest-frame UV. There is also uncer- 
tainty associated with neutral hydrogen absorption within the 
IGM. However, we model this on the average line-of-sight out to 
a given redshift, and our IGM absorption model will, therefore, 
not introduced any systematic effects. To investigate the effect 
that neutral host gas has on our results we re-ran our simulta- 
neous afterglow SED fits, but this time using only UV data red- 
ward of the Lyman forest (A > 1215A rest-frame). The best-fit 
results from these fits are summarised in Table [3] and the best- 
fit extinction curve is also plotted in Fig.[T](blue; long-dashed). 
The best-fit coefficients that describe the underlying UV contin- 
uum (ci and C2), and the height of the 2175A bump (C3), are not 
greatly affected by the removal of UV data blueward of 1215A, 
and remain consistent within their la errors. The new best-fit 
value of C4, which describes the far-UV extinction curvature 
at wavelengths A < 1700A, shows greater variation, increas- 
ing from 0.43^0;" in our full data-set analysis, to 1.79l|!);g4 
when data points that lie within the Lyman forest were removed. 
Nevertheless, the values remain consistent within 2a. 

Given that the removal of those data points blueward of 
A — I215A reduces the constraint on the far-UV curvature by 
over 50 it is not surprising that the C4 parameter shows 
the greatest variation in our fits with or without UV data at 
A < I215A. Furthermore, the fact that C4 increases in the 
latter fits suggests that neutral gas absorption within the host 
galaxy does not have a significant effect on our results. If the 
UV data blueward of 12 15 A had been significantly affected by 
host galaxy neutral hydrogen absorption, we would expect the 
omission of these data in our SED analysis to have yielded a 
smaller best-fit value of C4 relative to our initial analysis (i.e. the 
far-UV slope of the extinction law to have flattened), which is 
not the case. Furthermore, the best-fit extinction curve to our re- 
duced data set lies within the 3a confidence region of our best-fit 
results to the full data set (Fig. [T]i- This therefore suggests that 
our results are not significantly affected by additional absorption 
of the UV afterglow from host galaxy neutral hydrogen. 

5.2. Selection effects 

5.2.1 . How constrained is tine extinction curve slope? 

The optical slope of extinction curves is commonly param- 
eterised by the total-to-selective extinction parameter Ry— 
Ay I E{B — V), whereby for a fixed Ay, a larger Ry would im- 
ply a smaller amount of reddening, E{B — V), and thus a flatter 
extinction curve. As mentioned in section[3] our dataset does not 
have sufficient optical and NIR rest-frame coverage to constrain 
Ry. Nevertheless, our UV data do provide constraints on the 
value of Ry relative to the two coefficients, ci and C2, which de- 
scribe the continuum of the extinction curve at UV wavelengths. 
When we substitute ci for 2.09 — 2.84c2, the FM90 parameteri- 
sation that we have adopted (see Eq. [U then reduces to 

Ax^E{B-V)x (2) 
[Ry + 2.09 - 2.84C2 + C2/A + C3D{X-^) + aFiX'^)] 

for A in units of microns, and for A < 0.27^m, and where the 
other terms are as defined in section [3] It is clear from Eq. |2] 
that the parameters Ry and C2 both contribute to the wavelength 



^ Of the 46 data points at rest-frame A < 1700 A, where the C4 coef- 
ficient becomes relevant, 25 he within the Lyman forest (i.e. rest-frame 
A < I215A) 
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Fig. 3. 2-dimensional contourplot of C2 against Ry for C3 — 0.0 
(top) and C3 = 1.0 (bottom), and C4 = 0.4 in both cases. In 
both panels, the inner (black), middle (red) and outer (green) 
contours correspond to 68.3%, 90% and 99% confidence inter- 
vals, respectively. The best-fit ci and Ry values are indicated 
by a black cross, and the corresponding \Ry,c-2\ pairs for the 
mean SMC, LMC and Milky Way (MW) extinction curve are 
also shown. 



independent term. However Ry and C2 are not completely de- 
generate because the parameter also appears in the inverse 
wavelength dependent term. We have used XSPEC to determine 
what parts of Ry — parameter space can provide a statisti- 
cally acceptable fit to our GRB afterglow SED measurements. 
For this experiment we set the FM90 coefficients C3 and C4 to 
their best-fit values of 0.0 and 0.43 respectively. Both these pa- 
rameters also influence the UV slope of the extinction curve, 
and for now we therefore froze these in order to better illustrate 
the dependancy between Ry and c^. We explored the parameter 
range 2 < Ry < 6, 0.4 < C2 < 3.6, which covers the range 
in Ry th at has thus far b een measured (Gordon et al., 20031; 
Fitzpatri ck"& Massall2007h . The upper panel of Fig.[3]shows the 
best fitting [Ry, C2] pairs for the GRB SEDs together with their 
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Table 4. Best-fit results from FM90 model fits to the afterglow SEDs of 24 GRBs, all of which were all best-fit by a broken-power 
model in our initial SED analysis (see Table[T]l. 



GRB 


Local-dust" 








Power Law 




OLx ~ Oo 




model 


E{B - V) 


X^/dof 


Cl 


E(B - V) 


x'/dof 


(Y/N) 


050319 


mw 


^ ^ — 0.05 


58/78 


3.60+iii^^ 

'-'•^^—0.32 


^ — 0.10 


62/78 


Y 


050820A 


Imc 


0.06 ± 0.01 


177/143 


2.95 ± 0.25 


0.23 ±0.05 


190/143 


N 


060206 


mw 


< 0.08 


61/59 


4.39+?o';, 


< 0.03 


84/59 


N 


060418 


smc 


0.045+nE!7 

— U.UU I 


19/23 


— U.Ui 


< 0.01 


29/23 


N 


060729 


smc 


0.05 ±0.01 


176/191 




< 0.15 


180/191 


Y 


061121 


Imc 


0.17+^3 

— U.Uo 


87/93 


2.65+°-i° 

— U.4d 


0-73l°:li 


86/93 


N 


061126 


Imc 


— U.Uo 


170/140 


— U.4d 


0.36 ±0.05 


168/140 


N 


070110 


smc 


< 0.07 


129/112 


3.39lg;g^ 


0.82lO;l^ 


137/112 


Y 


080210 


Imc 


0.08 ±0.01 


23/25 


— l.oo 


0.09 ±0.02 


35/25 


N 


080411 


mw 


< 0.02 


81/83 


h 


0.00 


216/83 


Y 


080413B 


mw 




78/69 


2.94+g-33^ 


0.60lHo 


75/69 


Y 


080430 


smc 


< 0.07 


38/40 


o qi+0.86 


r, 91+0.09 

\j 1. Q Qg 


52/40 


Y 


080721 


mw 


0.13 ±0.05 


168/153 


o gi +0'31 


0.77 ±0.06 


168/153 


N 


080804 


Imc 


0.033 ±0.004 


38/42 


6 


< 0.01 


92/42 


N 


0809 16A 


smc 


0.21 ±0.05 


18/26 


1.24iro? 


0.64l°:l? 


19/26 


N 


081121 


smc 


n 023+° °"^ 

U.UZO_Q QQ4 


52/47 


6 


< 0.01 


92/47 


N 


081203A 


mw 


0.09 ±0.01 


76/80 


-0 76+'''^'' 


0.23 ±0.05 


61/80 


Y 


081222 


mw 


n 012+°""^ 

U.UiZ_Q QQ4 


34/46 


b 


< 0.01 


155/46 


N 


090102 


mw 


0.13 ±0.01 


22/27 


n 9O+0-66 
U-^^_0.83 


0.16 ±0.02 


35/27 


N 


090424 


mw 


0.37 ±0.03 


92/83 




0.52 ±0.08 


95/83 


Y 


090618 


smc 


p, f^ro+0.003 


24/18 


< 3.7 


< 0.02 


36/18 


Y 


091018 


smc 


0.033 ±0.005 


55/52 


b 


0.00 


101/52 


N 


091029 


Imc 


0.02 ±0.01 


23/23 


b 


< 0.02 


58/23 


Y 


091127 


Imc 


< 0.01 


186/173 


b 


0.00 


326/173 


N 



NOTE- All FM90 model parameters were fixed to the best-fit values from the simultaneous SED analysis (see TableO apart from the coefficient 
Cl, which determines the slope of the extinction curve. Column 8 indicates whether the optical and X-ray light curve decay rates are consistent 
(Y) or not (N) at epoch of the corresponding GRB afterglow SED. Those GRBs better fit by a single power-law spectral component and an FM90 
host extinction curve are highlighted in bold (see text for details). 

" The best-fit local dust model when using the broken power-law model, chosen from the mean SMC (smc), LMC (Imc) and Milky Way (mw) 

extinction laws. 

' Cl is unconstrained 



confidence contours. For comparison the average values of the 
[Rv, C2] parameters found along SMC, LMC and Milky Way 
sightlines are also shown. With C3 fixed at 0.0 we see that the 
mean {Ry, C2] values for the SMC and the LMC lie within the 
90% confidence region, whereas the Milky Way is outside of the 
99% confidence region. This finding is consistent with the dis- 
parity between the mean GRB extinction curve and mean Milky 
May extinction curve, as previously illustrated in the Fig. [T] 
curves. Fig. [3] confirms that for C3 — 0.0 and C4 = 0.43, we can 
rule out with 99% confidence that the mean GRB UV extinction 
curve is as flat as the mean Milky Way curve. 

To explore how much the presence of a 2 175 A bump in the 
GRB extinction curves would affect our constraints on Ry, C2, 
we repeated the previous experiment but with C3 set to 1.0 (with 
C4 fixed again to 0.43). The results are shown in the lower panel 
of Fig. [3] We see that, as before, the Milky Way [Ry, C2] param- 
eters he outside the outer confidence contour, and now the LMC 
[Rv, C2] parameters are also outside the 99% confidence region. 
Therefore, even if the mean GRB extinction curve does have a 
prominent 2175A bump, it is still incompatible at 99% confi- 
dence with an extinction law that is as flat as the mean Milky 
Way extinction law, and it is also incompatible with being as flat 
as the LMC extinction curve. 

Finally, we have also explored what constraints we can place 
on the Rv — C2 parameter space with C4 left as a free parame- 



ter, and we find that although with the extra free parameter the 
confidence regions widen, for both C3 — 0.0 and C3 = 1.0, the 
[Rv,C2] Milky Way values are still excluded at 99% confidence. 

We now explore further the slope of the GRB extinction 
curves at rest-frame UV wavelengths for the sample of GRBs 
that were excluded from our final sample. We use the parame- 
ter C2 (at a given Rv — 3.16) as our primary measure of the 
UV-slope, where A is in microns, and for A < 0.27/im. 

5.2.2. Are we biased against flat extinction curves? 

It is possible that our initial step of filtering out GRBs that appar- 
ently have a spectral break within their NIR to X-ray afterglow 
SEDs could preferentially reject sources with intrinsically flat 
dust extinction curves. This is because there is a chance that a 
flatter extinction law than that considered in our analysis may 
have been compensated by a best-fit optical spectral index that 
is flatter than the true value, as well as by a smaller total host 
reddening. 

To investigate if our method is biased against flat dust ex- 
tinction curves, we tested whether the GRBs with an apparent 
spectral break can alternatively be explained by unbroken power- 
laws with particularly flat dust extinction models. There were 24 
GRBs in our sample with an apparent spectral break, and we re- 
fit these with a single power-law continuum but using the more 
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flexible FM90 extinction curve prescription to model the host 
galaxy dust extinction. There are typically too few UV to NIR 
data points in these SEDs to constrain all the coefficients in the 
FM90 extinction model. The FM90 parameters were therefore 
fixed to the best-fit parameters from our simultaneous SED anal- 
ysis (Table [3]), with the exception of ci, which was left free to 
vary, and which is tied to C2. 

The results of these new fits are summarised in Table [H 
where we also give the goodness of fit to the afterglow SEDs 
provided by a broken power law model with normal (i.e. SMC, 
LMC or Milky Way) host galaxy dust extinction. A broken 
power-law with a local dust extinction model has the same num- 
ber of free parameters as a power-law with the FM90 extinc- 
tion curve, and so it is acceptible to compare directly the min- 
imum acheived by these two models. In only 5/24 cases 
(GRB 061121, GRB 061126, GRB 0804 13B, ORB 080721 and 
GRB 08 1203 A) a power-law with the FM90 extinction curve 
provided a better fit to the data than a broken power-law with 
dust extinction modelled on the local Universe. These five GRBs 
are highlighted in bold in Tableg] In the case of GRB 08 1203 A, 
the best-fit ci parameter lies within the best-fit value to the 
mean Milky Way and mean LMC extinction curve (see TableO. 
The improvement in the fit statistic relative to the SMC, LMC 
and Milky Way extinction law models is therefore likely to 
be due to the lack of a 2175A feature in our 'flat dust' ex- 
tinction curve model. Furthermore, we note that at the epochs 
where we have extracted SEDs for GRB 061121, GRB 061126 
and GRB 080721, the temporal decay rates measured in the 
X-ray band are different from those in the optical/NlR bands 
(GRB 061121; see fig. 10 of Page et al., 2007, GRB 061126; see 
fig. 4 & 7 of Gomboc et al., 2008 and see fig. 6 & 8 of Perley et 
al., 2008a, GRB 080721; see table 2 of Starling et al., 2009). The 
NIR through to X-ray afterglow emission for these three GRBs 
cannot therefore arise from the same emission component, and 
thus cannot lie on same NIR-X-ray spectral component. 

This therefore just leaves GRB 0804 13B'' as possibly having 
a host galaxy dust extinction curve that is flatter than those host 
extinction curves used in our SED analysis (section O. This 
GRB has a best-fit ci coefficient ci = 2.94lQ'3g, which is much 
larger than the best-fit FM90 coefficients corresponding to the 
mean SMC (ci = -4.47±0.19), LMC (ci = -1.28 ±0.34) and 
Milky Way (ci = 0.12 ± 0.11) values, and thus corresponds to 
a flatter extinction curve. This would suggest that only of order 
2% of GRBs that meet our initial parent sample selection criteria 
(49 GRBs; see section ITTT ) have flat dust extinction curves. 

A more conservative upper limit on the fraction of GRBs 
with a flat host extinction curve is given by considering the num- 
ber of GRB afterglows acceptably fit by a flat FM90 model, in- 
dependent of the goodness of fit of other dust models. For this, 
we only consider those GRBs with afterglows that decay at com- 
parable rates in the X-ray through to optical/NIR bands, for the 
same reason as was discussed above for GRBs 061 121, 061 126 
and 080721. I.e. it is only those GRB afterglows that have iden- 
tical temporal behaviour in the X-ray and NIR/optical/UV bands 
that can have a single NIR to X-ray spectral component. We 
therefore compare the temporal decay rates measured in the X- 
ray band and the NIR-optical bands for each of the GRBs in 



Tiie X-ray and optical/NIR afterglow light curves of GRB 0804 13B 
follow a complex evolution that do not satisfy t he closure relation s 
and are best explained by a two-component jet dFilgas et al.L 1201 Ih . 
However, within this model, the emission at the epoch of our SED is 
likely to be dominated by a single jet component, in which our single 
component spectral model would be valid. 



Table |4] to identify what fraction of these GRBs cannot have a 
single power law continuum. We find that 14/24 of these GRBs 
have decay rates in the X-ray band that differ from their decay 
rates in the NIR/optical bands, indicating that they must have a 
continuum that is more complex than a single powerlaw. This 
leaves 10 GRBs that can potentially have a single NIR to X- 
ray spectral slope. A column is included in Table |4] to indicate 
whether the optical and X-ray decay rates of each GRB are con- 
sistent at the epoch of the SED. Of these, 7/10 were accept- 
ably fitted by a power law and flat host galaxy extinction curve 
model at the 99% confidence level, although a further single 
GRB (GRB 090618) was significantly better fitted with a spec- 
tral break. We therefore consider a maximum of 6/49 GRBs, or 
an upper limit of 12% of GRBs in our sample to have flat host 
extinction curves. 

5.2.3. Redshift dependent extinction curves 

Our initial sample of 49 GRBs had a redshift range z = 0.5 — 4.0 
with a median redshift of (z) = 1.7, and in our final sample 
of 17 GRBs, the redshift range was from 0.7 to 3.1 with a me- 
dian of (z) =1.4. The median redshift of our GRB sample is 
lower than the median of (z) ^ 1.94 for the complete sample of 
Swift bursts^. This difference is the result of our selection crite- 
ria, which favours bright afterglows, and thus (relatively) low-z 
GRBs. It is currently not clear whether the galaxy dust extinc- 
tion curves change as a function of redshift for redshifts z < 5, 
and this bias in our sample should thus be kept in mind when 
applying the results presented in this paper to GRBs at redshifts 
z > 2. 

5.2.4. Host Ay dependent extinction curves 

A more obvious selection effect in our analysis is the bias against 
very dusty lines-of-sight. By requiring that the GRBs in our 
sample have both a spectroscopic redshift measurement and a 
UVOT u-band magnitude t; < 19 we remove GRBs from our 
sample with highly extinguished afterglows, both due to dusty 
regions within their host galaxy and/or high Galactic extinc- 
tion along their line-of-sight. If the extinction law of GRB host 
galaxies varies with the host galaxy dust column density, then 
our analysis will clearly only be applicable to GRBs with a host 
galaxy visual extinction smaller than a certain threshold, above 
which the GRB afterglow would be too extinguished to enter 
our sample. There is empirical evidence that the prominence of 
the 2 175 A extinction bump is related to the total visual extinc- 
tion. Ay, along the line-of-sight. The host extinction curves of 
GRB 070802, GRB 080605, GRB 080607 and GRB 080805 are 
the only ones with a spectroscopically detected 2175A extinc- 
tion bump, and the afterglows are also of the most heavily extin- 
guished. This point is also discussed in lZafar et al.l (1201 Ih . 

Such a relation between the prominence of the 2 1 75 A feature 
and visual extinction. Ay, is not necessarily indicative of evolu- 
tion in the dust extinction curve, and may be a result of there 
not being sufficient extinction in low dust environments for the 
2175A feature to be detected with any statistical significance. It 
is, nevertheless, possible that the results presented in this paper 
are only valid for the host galaxies of those GRBs that have a 
total visual extinction Ay < 1. In view of this prospect, in the 
next section we investigate further the origin of differences in 
the host galaxy extinction curves of GRBs with only moderately 
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extinguished afterglows, such as in our sample, and those that 
are highly extinguished (i.e. Ay > 1). 

5.3. Variation in GRB liost extinction laws 

The improved positional accuracy of GRBs provided by the 
rapid response of Swift and (semi-)robotic ground-based tele- 
copes, in particular those equipped with NIR instruments, has 
significantly improved our capabilities over the last half decade 
to study highly dust-extinguished GRBs and their host galax- 
ies. Dedicated host galaxy follow-up programmes of highly 
dust-extinguished GRB candidates are revealing a sample of 
GRB host galaxies that are typically more massive, luminous, 
and chemically evolved than the typical host galaxies of rela- 
tively unextinguished GRBs (GRB 020127; Berger et al. 2007, 
GRB 030115; Levan et al. 2006, GRB 070802; Ki-uhler et al. 
201 1, GRB 080325; Hashimoto et al. 2010, GRB 080607; Chen 

etal. 2011). 

In the case of GRB 080607. I^ochaska et al.1 (|2009|) argued 
the majority of the afterglow dust extinction to come from dust 
within an intervening molecular cloud. However, it is notable 
that the characteristic properties of the host galaxy a re also atyp- 
ical w hen compared to optically- selected samples (Chen et al.l, 
1201 lb . The host galaxy of GRB 080607 had a stellar mass 
M^, ^ 8 X 10^ Mq, which is almost an order of magnitude 
larger than the mea n stellar mass o f optically selected samples 
((M*) ~ 10^ M^: ISavaglioltaLL 12009; Kruhleret al., 2011), 
and it was very red, with R — K > 5 (iChen et al.LI201 1) . This 
is also the only GRB to have a robust detect ion of molecular 
hydro gen absorption in its afterglow spectrum (Prochaska et al.L 
I2009h . 

On the other hand, GRB 070306 and GRB 100621 A were 
heavily extinguished GRBs (Ay ^ 5.5 and Ay ^ 3.8 re- 
spectively), but both had very blue host galaxies, with R — K 
colours compar able to the ho s t galaxies of relatively une xtin- 
guished GRBs (iJaunsen et all 120081; iKruhler et all l201lh . in- 
dicative of a very clumpy distr ibution of dust. In a sample of 
six highly extinguished GRBs, Perle v et"an (12009') also found 
evidence for a patchy dust distribution. Furthermore, although 
GRB 070306 had one of the largest ste llar masses measure d for 
a long GRB host galaxy (2 x 10^° M^; lKrUhler et all 1201 ih . the 
host stellar mass for GRB 100621 (10^ Mq) was comparable 
to that of optically bright afterglow host galaxies. The relation 
between afterglow dust extinction and host galaxy properties, 
including the host extinction curve, is therefore not clear, and 
further investigation of the host galaxies of heavily extinguished 
afterglows will be needed to address these issues. 

6. Conclusions 

Using a sample of 17 GRBs with good wavelength coverage 
and a single NIR through to X-ray spectral emission component, 
we have been able to constrain the host galaxy dust extinction 
properties as a function of wavelength for GRBs with total host 
galaxy visual extinction Ay < 1. By studying the host galaxy 
dust extinction properties for a sample of GRBs rather than for 
a single, well-observed, case, the analysis presented in this pa- 
per allows us to investigate the average dust extinction proper- 
ties in the host galaxies of GRBs, albeit only for those with a 
moderate host galaxy visual extinction (Ay < 1.0). Using the 
[pitzpatrick & Massa (1990) analytic extinction curve model to 
fit our afterglow sample, we find the best-fit to have a slope 
blueward of 2000A that is intermediate to that of the SMC 



and LMC mean extinction curves, and it has no 2 175 A extinc- 
tion bump. Within the uncertainties, the best-fit extinction law is 
consistent with both the SMC and LMC mean extinction curves. 
However, our data are inconsistent with having a host galaxy ex- 
tinction curve as flat as the average Milky Way curve blueward 
of - 2000A. 

Our best-fit FM90 parameterisation of the GRB host extinc- 
tion curve provides an acceptable fit (at 3a confidence) to a full 
sample of 49 afterglow SEDs that we initially analysed, provid- 
ing typically a similar or improved fit to the afterglow SEDs 
compared to fits that use local dust extinction curve models. Of 
those GRBs with host galaxy dust extinction measured at 90% 
confidence, around 10% were, however, better fit by an extinc- 
tion curve that contained the 2 175 A dust feature. Furthermore, 
up to 12% of the 49 afterglow SEDs analysed could have a host 
galaxy extinction curve that is significantly flatter at rest-frame 
UV wavelengths than the mean Milky Way extinction curve. For 
this 12% of GRBs, the improved fit provided by a flatter UV- 
extinction curve model increases the host galaxy total extinc- 
tion relative to fits with local extinction curve models. This is 
compatible with the trend increasingly observed between sig- 
nificantly extinguished GRBs (A„ ^ 1) and flat host galaxy 
extinction curves. 

The relatively UV-steep host galaxy extinction curve for our 
sample of modestly extinguished GRBs, and the much flatter ex- 
tinction curves with more pronounced 2 175 A features of host 
galaxies of more extinguished GRBs, is adding to the empiri- 
cal evidence that the extinction curve of GRB hosts vary as a 
function of host dust abundance and/or some other galaxy prop- 
erty. In view of this, a future test would be to study the aver- 
age GRB host extinction curve properties for a sample of GRB 
afterglows with litde (e.g. Ay < 0.3) and high (Ay > 1.0) 
host galaxy extinction. Such analysis would require a sample of 
bright, relatively unextinguished GRBs in the former case, and 
a sample of extinguished GRBs with good broadband rest-frame 
NIR through to X-ray coverage in the latter, and preferably with 
no spectral break within the observed bandpass. The current size 
of such a sample remains relatively small. However, examples of 
highly extinguished cases such as GRB 080607 show that such 
GRB samples are attainable with time, as long as rapid-response 
observing programmes that cover the optical through to NIR af- 
terglow remain available. 
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Fig.A.l. NIR through to X-ray SEDs of the 17 GRBs in our final sample with the best-fit FM90 model overplotted. Frequency 
is plotted in the rest-frame and flux density is in arbitrary units for clarity. The SEDs have been corrected for Galactic reddening 
and soft X-ray absorption, and both the host galaxy absorbed/extinguished (solid) and intrinsic (dashed) best-fit afterglow SED are 
indicated. The afterglow SEDs are ordered by afterglow host galaxy extinction, with the most extinguished at the bottom, and the 
least at the top. 



18 



P. Schady et al.: The Dust Extinction Curve of Gamma-Ray Burst Hosts 



16 
14 
12 



(/3 

H 

2 8 



6 - 



4 - 



2 - 



- 



-2 - 



-i — I I I I 1 1 1 1 1 — I I I I 1 1 1 1 1 — I I I I 1 1 1 

080210 Av=0.26 ^^fez - 
(bknp/LMC) ^fer 

08 1 203 A Av=0.29 fes: - 
(bknp/MW) 

050319 Av=0.32 
rbkno/MW) 

080721 Av=0.39 
(bknp/MW) 



090102 Av=0.41 
(bknp/MW) 

061121 Av=0.55 
(bknp/LMC) 

061126 Av=0.57 
(bknp/LMC) 

0809 16A Av=0.62 
(bknp/SMC) 



-| 1 1 I I I I I I 1 1 1 I I I I 1 1 1 1 1 I I I I 1 1 




08041 3B Av=0.84 
(bknp/MW) 



090424 Av= 1.1 4 
(bknp/MW) 





_l I I I I I I I I I I 



10 



14 



J5 



J6 



M 



10'" 10'" 10" 10 
Restframe Frequency (Hz) 



18 



10 



19 



Fig.A.2. NIR through to X-ray SEDs of the 32 GRBs analysed that did not enter our final sample. Frequency is plotted in the 
rest-frame and flux density is in arbitrary units for clarity. The SEDs have been corrected for Galactic reddening and soft X-ray 
absorption, and both the host galaxy absorbed/extinguished (solid) and intrinsic (dashed) best-fit afterglow SED are indicated. The 
afterglow SEDs are ordered by afterglow host galaxy extinction, with the most extinguished at the bottom, and the least at the top. 
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Fig. A.2. (continued) 
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Fig. A.2. (continued) 



